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We discuss a search for CPT violation in neutral charm meson oscillations. The 
data come from the Fermilab fixed-target experiment FOCUS. While flavor mixing 
in the charm sector is predicted to be small by the standard model, it is still possible 
to investigate CPT violation through study of the proper time dependence of the 
asymmetry in right-sign decay rates for D° and D°. Using present limits for 
D° — D° mixing we infer bounds on charm CPT violation using data from FOCUS. 


1 Introduction 

The combined symmetry of charge conjugation (C), parity (P), and time re¬ 
versal (T) is believed to be respected by all local, point-like, Lorentz covariant 
field theories, such as the standard model. However, extensions to the standard 
model based on string theories do not necessarily require CPT invariance, and 
observable effects at low-energies may, be within in reach of experiments study¬ 
ing flavor oscillationsEl A frameworlfl in which indirect CPT and T violating 
parameters appear has been developed which allows experimental investigation 
in neutral meson systems. 

For quite some time searches for CPT violation have been made in the neu¬ 
tral kaon system in which particle-antiparticle mixing is large. KTe\Q reports a 
bound on the CPT figure of merit t'k = \rriK° — m^o\/m^o < (4.5±3) x 10 -19 . 
CPT tests in B° meson decay have been made by OPAIb at LEP, and by Belle 
at KEK which has recently reported! tb = \rriB0 — < 1-6 x 10 -14 . 

To date, no experimental search for CPT violation has been made in the 
charm quark sector. This is probably due in part to the expected suppression 
of D° — D° oscillations in the standard model, and the lack of a strong mixing 
signal in the experimental data. Recent mixing searches include a study of life¬ 
time differences between charge-parity (CP) eigenstates from FOCUS, which 
reported! a value for the parameter yep = (3.42 ± 1.39 ± 0.74)%. The CLEO 
Collaboration has reported 95% confidence level bounds on mixing parameters 
x' and y ' (related to the usual parameters x and y by a strong phase shift) 0 
(l/2)a/ 2 < 0.041% and —5.8% < y’ < 1%. FOCUS has reported! a study of the 
doubly Cabibbo suppressed ratio ( Rdcs ) f° r the decay D° —► and has 
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extracted a contour limit on y' (of order few %) under varying assumptions of 
Rdcs and x’ . The question arises - what can be learned about indirect CPT 
violation given the apparent smallness of mixing in the charm system? It turns 
out that even in the absence of a strong mixing signal one can still infer the 
level of CPT violation sensitivity through study of the time dependence of D° 
decays, which we show in this paper. 


2 Proper Time Asymmetry 

Time dependent decay probabilities into right-sign and wrong sign decay modes 
for neutral mesons which express the CPT violation have been developed in 
a general frameworkEl For the decay of D° to a right-sign final state / (which 
could be a semileptonic mode, or a Cabibbo favored hadronic mode), the time 
dependent decay probability is 

p f (t) = \{f\T\D°m 3 

= -|F| 2 exp(—^t) x [(1 + |£| 2 )coshA 7 t /2 + (1 — |£| 2 )cosAmt 
—2 Re £ sinhA 7 f /2 — 2 Im £ sin Amt]. (1) 

The time dependent probability for the decay of D° to a right-sign final state 
/, Pf{t), may be obtained by replacing in the above equation £ —► — £ and 
F —> F. In the formula, F represents the basic transition amplitude for the 
decay D° —> /, A 7 and Am are the differences in physical decay widths and 
masses for the propagating eigenstates and can be related to the usual mixing 
parameters x = A M/T and y = AT/2T. The complex parameter £ controls 
the CPT violation and is seen to modify the shape of the time dependent decay 
probabilities. Expressions for wrong-sign decay probabilities involve both CPT 
and T violation parameters which scale the probabilities, leaving the shape 
unchanged. Using only right-sign decay modes, the following asymmetry can 
be formed, 


Acprit) 


Pf(t) - P(t ) 
p f -(t) + p(ty 


which is sensitive to CPT violating parameter £: 


Acpt 


2 Re £ sinhA7t/2 + 2 Im £ sin Amt 
(1 + |£| 2 )coshA 7 f /2 + (1 — |£| 2 )cosAmt 


( 2 ) 

( 3 ) 


We can gain insight into the anticipated experimental sensitivity by plot¬ 
ting these functions with some reasonable assumptions. We use 95% confidence 
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level upper bounds on the mixing parameters x and y of 5%, which is at the 
upper range of the current experimental sensitivity, as discussed previously. In 
Fig. 0(a) we plot the proper time decay probabilities for D° decay under the 
assumption of CPT violation at the level of Re £ =5%, Irn £ = 5%, which are 
independent parameters in the framework. One sees a CPT violation-induced 
wrong-sign contribution which vanishes at zero proper time and at long proper 
times. This causes a distortion from a purely exponential decay of a D° (and 
D°), which is then visible in the asymmetry plot, Aqpt as shown in Fig. 0(b). 
Because of the small oscillation frequency and short lifetime, one sees only the 
start of the oscillation, growing beyond 0.3% at long proper times. Evident 
from Eqn. 0 is that positive values of the Re £ and Im £ work to oppose one 
another in the asymmetry in a linear fashion. This is shown in the nearly 
linear behavior of Acpt in Figs. 0(c,d) which have Im £ = 5%, Re £ = 0 
and Im £ = 0, Re £ = 5% respectively, and consequently CPT asymmetries 
larger by a factor of 10 at long proper times . In practice, experiments will be 
sensitive to either Re £ or Im £, but not both simultaneously. 


3 Data Analysis 

In this paper we investigate the current experimental sensitivity for a CPT 
violating signal using data collected by the FOCUS Collaboration during the 
1996-97 fixed-target run at Fermilab. FOCUS is an upgraded version of the 
E687 spectrometer. Charm particles are produced by the interaction of high 
energy photons (average energy w 180 GeV for triggered charm states) with 
a segmented BeO target. In the target region charged particles are tracked 
by 16 layers of microstrip detectors. These detectors provide excellent vertex 
resolution. Charged particles are further tracked by a system of five multi¬ 
wire proportional chambers and are momentum analyzed by two oppositely 
polarized large aperture dipole magnets. Particle identification is determined 
by three multi-cell threshold Cerenkov detectors, electromagnetic calorimeters, 
and muon counters. 

The data analysis is as follows. We analyze the two right-sign hadronic 
decays D° —> K~n + and D° — > K + n~. Contamination by doubly Cabibbo 
suppressed (wrong-sign) decay occurs but is small owing to the small branching 
ratio, and its effect on the right-sign signal yield estimates in this analysis is 
negligible. We use the soft pion from the decay D* + —> D°tt + to tag the flavor 
of the D at production, and the kaon charge in the decay D° —► K~tt + to 
tag the D flavor at decay. (Charge conjugates are assumed in this paper.) 
D° —+ K~-k + events were selected by requiring a minimum detachment of 
the secondary (decay) vertex from the primary (production) vertex of 5 at- 
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Figure 1: (a) Proper time decay probabilities for D° decay, right-sign (solid) and wrong-sign 

(dashed) and (b) Acpt with Re £ = 5%, Im £ = 5%. (c) Acpt with Re £ = 0, Im £ = 5%. 
(d) Acpt with Re £ = 5%, Im £ = 0. 


The primary vertex was found using a candidate driven vertex finder which 
nucleated tracks about a “seed” track constructed using the secondary vertex 
and the reconstructed D momentum vector. Both primary and secondary 
vertices were required to have confidence level fits of greater than 1%. The 
Z?*-tag is accomplished by requiring-.the D* — D° mass difference to be less 
than 3 MeV/c 2 of the nominal valueu 

Kaons and pions were identified using the Cerenkov particle identification 
cuts. These cuts are based on likelihood ratios between the various stable 
particle hypotheses, and are computed for a given track from the observed 
firing response (on or off) of all cells within the track’s (/3 = 1) Cerenkov light 
cone in each of three multi-cell, threshold Cerenkov counters. The probability 
that a given track wil fire a given cell is computed using Poisson statistics based 
on the predicted number of photoelectrons striking the cell’s phototube under 
each particle identification hypothesis; the accidental firing rate of the cells, 
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which depended on beam intensity, was also used. The product of all firing 
probabilities for all cells within the three Cerenkov cones produces a x 2 -like 
variable called W t = - 2xlog(likelihood) where i ranges over electron, pion, 
kaon and proton hypotheses. For the K and the ir candidates we require Wi to 
be no more than 4 greater than the smallest of the other three hypotheses ( Wi — 
Wmin < 4) which eliminated candidates that highly likely to be misidentihed. 
In addition, D° daughters must satisfy the slightly stronger Kir separation 
criteria W n — Wk > 0.5 for the K and Wk — W n > —2 for the n. Doubly 
misidentihed D° —> K~ "7r + is removed by imposing a hard Cerenkov cut on the 
sum of the two separations {{W^—W^k+^Wk—W^)^ > 8) when the invariant 
mass of the track pair with the K and tt hypothesis exchanged is within 4a 
of the nominal D° mass. Fig. || shows the invariant mass distribution for two 
D*-tagged, right-sign decays D° —> Iv - 7r + and D° —► K + ir~. Approximately 
49700 signal events were used in the analysis described below. 



Figure 2: Invariant mass distribution for the sum of D° and D° right-sign decay candidates. 

The reduced proper time is a traditional lifetime variable used in fixed- 
target experiments which use the detachment between the primary and sec¬ 
ondary vertex as their principal tool in reducing non-charm background. The 
reduced proper time is defined by t' = (£ — Na^/^^c) where i is the distance 
between the primary and secondary vertex, at is the resolution on l, and N is 
the minimum detachment cut required to tag the charmed particle through its 
lifetime. Fig. shows reduced proper time distributions for the two right-sign 
decays D° —> I\~tt + and 5° —> K + ir~. 

We plot the difference in right-sign events between D° and D° in bins of re¬ 
duced proper time t' in Fig. |J. For each data point, the background subtracted 
yields of right-sign D° and D° are used in forming the ratio. In the absence 
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Figure 3: Background subtracted reduced proper time distributions for D° and D°. 


of detector acceptance corrections, this is equivalent to Acpt as defined in 
Eqn. |^. Monte Carlo studies showed that detector and resolution effects cause 
less than 5% changes in the t' distribution as measured by deviations from a 
pure exponential decay. Acpt is additionally insensitive to detector effects 
since they tend to cancel in the ratio. 

Because of the QCD production mechanism for photoproduced charm 
mesons, more D° than D° are produced in the FOCUS data sample. This has 
been previously investigated in photoproduction by E687, in which the produce 
tion asymmetries were studied in the context of a string fragmentation modellij 
The effect in the Acpt distribution is to add a constant, production-related 
offset, which can be accounted for in fits to the asymmetry distributions. 

We infer the CPT sensitivity for a fixed set of experiment bounds on 
mixing without CPT violation. The Acpt data in Fig. are fit using the time 
dependent probability functions for D° and D°, as defined by Eq. [l] for the case 
of D° decay. The allowed fit parameters are a constant production asymmetry 
parameter and Re £ (fits with Im £ gave similar results for sensitivity). In the 
fit shown we used values for the x and y mixing parameters of 5%, and infer 
one standard deviation errors on either Re £ or Im £ of approximately 10%, 
corresponding to 95% confidence level upper bounds of roughly 25%. 

In the CPT and Lorentz violating extension to the standard model@ the 
CPT violating parameters may depend on lab momentum, orientation, and 
sidereal timeHa They are expressed in terms of Lorentz violating coupling 
coefficients which depend on the flavor of the valence quark states. (For this 
reason, CPT violation in the K , D, and B systems need not be same.) In the 
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Figure 4: Asymmetry as a function of reduced proper time (left). At the right is plotted the 
variation of the fitted CPT parameter Re £ as a function of local sidereal hour. 


case of FOCUS, a forward, fixed-target spectrometer, the £ parameter assumes 
the following form 

0,P) = ^^[Aa 0 +/3Aa z cosx 

+/3sinx(AaysinOt + AaxCOsUt)], (4) 

where SI and t are the sidereal frequency and time respectively, and X , Y, Z are 
non-rotating coordinates with Z aligned with the Earth’s rotation axis. Our 
sensitivity quoted above should then be regarded as that for |£| averaged over 
momentum and sidereal time. From this quantity one can infer the sensitivity 
for the Lorentz violating charm and u quark coefficients: Ado + 0.6 x Adz < 
10 -15 GeV, where we have used 53 degrees as the angle between the FOCUS 
spectrometer axis and the Earth’s pole. 

We searched for a sidereal dependence, which would be sensitive to co¬ 
efficients Aa.v and A ay, by dividing our data sample into four-hour bins in 
sidereal time and repeating our fits. The resulting distribution, shown also in 
Fig. U, can be fit using Eq. [| above to yield bounds on Aa.v and A ay of similar 
sensitivity. 

4 Conclusions 

We have investigated the sensitivity to a CPT violating signal in neutral charm 
meson oscillations, the first search of its kind in the charm system. Our prelim¬ 
inary study using FOCUS data indicates 95% confidence level upper bounds 
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on charm CPT violation parameters Re £ and Im £ of approximately 25%, 
corresponding to bounds on the CPT and Lorentz violating charm coefficients 
of < 10” 15 GeV. 
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